Kinematic reduction of reaction-diffusion fronts with multiplicative noise: 
Derivation of stochastic sharp-interface equations 
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We study the dynamics of generic reaction-diffusion fronts, including pulses and chemical waves, 
in the presence of multiplicative noise. We discuss the connection between the reaction-diffusion 
00 . Lan gevin-like field equations and the kinematic (eikonal) description in terms of a stochastic moving- 

boundary or sharp-interface approximation. We find that the effective noise is additive and we relate 
(— I I its strength to the noise parameters in the original field equations, to first order in noise strength, 

. but including a partial resummation to all orders which captures the singular dependence on the 

microscopic cutoff associated to the spatial correlation of the noise. This dependence is essential for 
a quantitative and qualitative understanding of fluctuating fronts, affecting both scaling properties 
and nonunivcrsal quantities. Our results predict phenomena such as the shift of the transition point 
between the pushed and pulled regimes of front propagation, in terms of the noise parameters, and 
, the corresponding transition to a non-KPZ universality class. We assess the quantitative validity 

of the results in several examples including equilibrium fluctuations, kinetic roughening, and the 
noise-induced pushed-pulled transition, which is predicted and observed for the first time. The 

S analytical predictions are successfully tested against rigorous results and show excellent agreement 
with numerical simulations of reaction-diffusion field equations with multiplicative noise. 

£5 . PACS number(s): 05.40.-a,05.45.-a,82.40.Ck,47.54.+r 

O ' 

I. INTRODUCTION 

>' 

The dynamics of localized solutions in the form of fronts or pulses in reaction-diffusion systems has received a great 
deal of attention for a long time in the context of nonequilibrium extended systems |lj . Examples of fronts formed 
by stable regions invading unstable or metastable ones are found in a large variety of physical, chemical or biological 
systems, and have been studied in great detail [2-[|. In the context of excitable media and chemical waves, extended 
pulses do also exhibit a rich phenomenology jf^-E]. As opposed to fronts, excitation waves or pulses are such that the 
region behind them eventually returns to the same linearly stable state that is ahead. In this case, in dimensions 
higher than one there may be open ends of the pulse region which give rise to spiral (2d) or scroll (3d) waves. There 
has been much interest in the study of such objects from the fundamental point of view of pattern forming dynamics, 
but also because of their potential applications in biological systems, such as in cardiac tissue JTo) ] . In this paper we 
will not consider the case of open ends, so unless otherwise specified, we will refer indistinctly to fronts and pulses 
P^j ' under the common term of 'fronts'. 

One aspect which has received increasing interest in recent years has been the effect of fluctuations both of internal 
and external origin in the dynamics and in the roughening properties of fronts p]-^9[. More recently the effect 
^ , of noise in pattern forming dynamics of chemical waves has been fostered by the development of the experimental 
capability to introduce external spatio-temporal noise in a controlled way in different photosensitive nonlinear chemical 
reactions, through the optical projection of computer-designed spatio-temporal fluctuations in the local illumination 
conditions, acting as a stochastic control parameter j2^23|. From a theoretical point of view, a common starting 
point to study fluctuations is in terms of master equations defining evolution of reacting and diffusing particles in 
a lattice pd| ]. The connection of this microscopic level of description to the mesoscopic one in terms of Langevin 
field equations has proven a rather subtle issue, in connection with the distinction of the so-called pushed vs pulled 
fronts Only recently a complete understanding of the instrinsically different nature of these two types of fronts, 
and the correspondingconsequences concerning the presence of cutoffs [^J-^, and the effects of fluctuations has 
started to emerge [ [L 2| , pL 3|Jl 5| — |l 7|| . It has been shown, for instance, that pulled fronts define a new universality class of 
kinetic roughnening different from the Kardar-Parisi-Zhang universality class Jp^,[^]. On the other hand, it has been 
shown that intrinsic noise at the microscopic level may induce a morphological instability at the macroscopic level of 
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description |TJ]. In this paper we will be mostly concerned with the macroscopic description of pushed fronts with 
fluctuations, but also on how this description incorporates the transition to pulled fronts induced by the noise itself. 

In the absence of noise, and in the appropriate limit, the description of fronts and pulses defined by reaction-diffusion 
field equations can be reduced by means of a moving boundary approximation to a kinematic description in terms of 
much simpler local equations |9|p8|]. This procedure, which is mathematically well grounded in the framework of the 
so-called inertial manifold reduction, can be carried out systematically, for instance, through asymptotic matching 
techniques using the front thickness as a small parameter defining a singular perturbation problem p9[. This is 
commonly used to relate macroscopic interface equations to order-parameter or phase-field descriptions in many 
problems involving interface dynamics, such as in solidification pO 31|, viscous fingering p2] etc. In the case of 
chemical pulses in excitable media in the limit of weak excitability, this leads to a local equation where the normal 
velocity of the pulse is a constant plus a correction proportional to curvature ||. This is often called the eikonal 
equation. In the case of open ends, a similar local equation can be derived for the motion of the endpoint of the pulse 
f^H- For pushed fronts it can be shown that for smooth, long wavelength deformations of the front, the separation 
of time scales between the soft deformation modes and the internal degrees of freedom of the fields leads naturally 
to the same eikonal equation. For pulled fronts, however, this separation of time scales does not exist, the relaxation 
being algebraic instead of exponential, and a local moving-boundary approximation is not justified [p9f . 

The kinematic description in terms of eikonal- like equations is a very useful approximation from both a theoretical 
and a practical point of view. In the context of the study of universality of fluctuation properties p3| , |34| for instance, it 
leads naturally to the relevant effective universal description of a broad class of systems in terms of the KPZ equation 
|j35f . It is also very useful for numerical simulation purposes, since it avoids resolving the fine structure of the front 
and the bulk degrees of freedom, which become irrelevant, dealing only with the kinematic degrees of freedom of an 
object of reduced dimension. 

When noise is present in the original field equations, however, the situation is not so clear. Stochastic eikonal 
equations have proven useful in a phenomenological description of the dynamics of pulses and spiral waves in pho- 
tosensitive chemical systems with external noise imposed on the illumination conditions pl| , p2[ . Such description, 
however, relied on some fitting parameters and some uncontrolled hypothesis on the way the noise must enter the 
kinematic equations. These results, together with the fact that the statistical properties of the noise present in those 
experiments are fully controlled, clearly call for a more 'microscopic' derivation of stochastic kinematic equations 
corresponding to Langevin reaction-diffusion field equations with multiplicative noise, with no free parameters. The 
connection between bulk and interface fluctuations has been worked out so far only for equilibrium fluctuations in 
coarse-grained, Ginzburg-Landau-like equations 36-3^]. In such cases, the sharp- interface limit can be taken at the 
level of the free-energy, and the existence of a fluctuation-dissipation theorem then yields the proper way to incorpo- 
rate thermal fluctuations into the effective interface equations. However, in many nonequilibrium systems, for instance 
in the context of reaction-diffucion problems, a free-energy or generically a Liapunov functional may not exist and 
no fluctuation-dissipation relation may be invoked for external noise. In such cases the connection between the bulk 
description and the effective interface description in presence of fluctuations must be worked out at the level of the 
dynamical equations. 

The purpose of this paper is to address this point by deriving stochastic eikonal equations, including the complete 
specification of the noise statistics in connection to that of the noise in the mesoscopic field equations. In particular 
we shall focus on the singular dependence on the spatial cutoff when noise is multiplicative, and its importance in the 
quantitative description of the statistical properties of the front fluctuations. The predictions will be tested against 
numerical simulations of reaction-diffusion equations, and also in cases where exact results are available concerning 
the spectrum of interface fluctuations. We will also see that the stochastic eikonal equation derived is consistent with 
the scenario of the pushed-pulled transition, and that, changing the spatial cutoff or the noise intensity may have 
effects such drastic as changing the universality class of kinetic roughening of the front through that transition. 

Although our obtention of the stochastic eikonal equation here presented is not a first-principles rigorous derivation, 
we will provide sufficient evidence to conclude that the result is the correct one to lowest order in noise intensity, 
including the singular dependence on the noise correlation cutoff (which involves a partial resumation to all orders), 
within the long time and length scales limits which are inherent to the kinematic descrition itself. However, it does 
not apply to situations in which the front dynamics is nonlocal, such as solidification fronts or viscous fingers, where 
a different type of formulation is appropriate even in the deterministic case PQ-p3,E0[ . 
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II. KINEMATIC REDUCTION FOR GENERIC REACTION-DIFFUSION SYSTEMS 



Let us consider a vectorial field </>(x, t) with N components </>(x) = </>i(x), 4>n(x.) in a d— dimensional space with 
x = X\, Xd, which obeys a reaction-diffusion equation with multiplicative noise of the form: 

^=i3v 2 + f(0)+ £ 1 / 2 g (0)r / (x,t), (1) 

where rj(x.,t) is a gaussian noise with zero mean and correlation given by 

(77(x, t) v ( x ', 0) = 2A- d C(|x - X '\/X)6(t - t'). (2) 

We take a one-component noise for simplicity, as the natural case when it originates in fluctuations of a single control 
parameter. The generalization of the formalism to multicomponent noise is straightforward. Notice the asymmetry 
with which we treat the spatial and temporal correlator of the noise. As we will see, this reflects a nontrivial issue 
related to the intrinsically different nature of the white noise limit in space as opposed to time. We have taken in Eq. 
(^) the noise as delta-correlated in time. This temporal white noise limit is well behaved, once a prescription for the 
multiplicative noise term in Eq. (|l|) has been chosen. For external fluctuations the physically relevant prescription is 
to consider the white noise as the limit of some properly defined correlated noise. This corresponds to the well-known 
Stratonovich prescription pi] . On the other hand, the spatial noise must always be defined as colored, its white 
(uncorrelated) limit being intrinsically ill-defined. Hence the notation with the function X~ d C(r / A) in the correlator, 
meaning a general correlation function, dependent on some correlation length A, which in the limit A — > is such that 
\~ d C{r / X) — > 8(r). The fact that the spatial continuum limit A — > does not exist in the Stratonovich interpretation 
is thus reflected in the fact that, even if A is much smaller than any other length scale in the problem, the existence 
of such microscopic cutoff always shows up in the quantitative description of the large scale behavior and cannot 
be reabsorbed in a redefinition of parameters. That is, the noise cannot be considered as effectively white in space, 
regardless of how small the noise correlation length is. 

The multiplicative noise term present in Eq. (|l|) has an average value different from zero. Applying Novikov 
theorem jl2) , we get 

e^g^Mx,*)) =eA- d C(O)(G(0)) (3) 

where 

G ^)-E%r^(0). (4) 

This suggests to separate the average contribution from the multiplicative noise term and rewrite Eq. ([!]) in terms of 
a renormalizcd potential and a zero average noise, 

^=^V 2 + h(0)+ £ 1 /2 O(0;Xji)j (5) 

where 

h(0) =f(0)+eA- d C(O)G(</>) (6) 

and 

x, t) = G(0)j7(x, t) - ei A- d C(O)G(0). (7) 

Here the new noise ft has zero average. Notice that this decomposition corresponds to transform Eq. ([!]) into its 
equivalent Ito stochastic equation in the white noise limit. Accordingly, the stochastic term f2 reduces to 

lim n{<f>, x, t) = G(0)t7/(x, t) (8) 

A — >0 

where r\i (x, t) is now a white noise in the Ito interpretation. The deterministic term h includes thus noise effects 
through the so-called Stratonovich term that has been added to f . These noise effects on the deterministic part of 
Eq. (||) depend on a 'dressed' noise intensity E\ which contains the singular dependence on the the spatial cutoff in 
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the form of the 'bare' e as £\ = eC(0)A~ d . The important point is that the spatial white noise limit in the continuum 
equations is mathematically well-defined for an Ito noise. This has been proven rigorously for relatively broad classes 
of equations (see for instance Refs. On the contrary, this cannot be the case for a Stratonovich noise, as it 

is obvious from the singular dependence on A. The practical implications of these fact are that, while the singular 
contribution of A in the Stratonovich term must be kept explicitly, the dependence on A contained in £1 will indeed 
be weak (nonsingular) , and indeed negligible if A is much smaller than the other length scales in the problem, in 
particular the front thickness. Therefore, the splitting of the Stratonovich noise by means of a deterministic term plus 
an Ito noise has the virtue to isolate the singular dependence on the microscopic cutoff. This will also define a useful 
partial resummation of orders in ex- After this decomposition, the theory will thus be correct to first order in e, to 
all orders in e/X and to order A for the regular dependence of ft in the noise correlation length. 

Eq. (g) will be our starting point. We are interested in representing the dynamics of a front obeying Eq. rt5ft as 
the evolution of a d — 1 surface. In this effective dynamics we assume the details of the front structure (at scales of 
order of or smaller than the front thickness) to be unimportant. In what follows we will consider the evolution of a 
1-dimensional front embedded in a 2-dimensional system. In the procedure we are going to apply we will write the 
evolution equation @ in the curvilinear coordinates defined by the Id front in the sharp-interface limit, and obtain 
the evolution equation for this front as a solvability problem with the basic assumption that curvature and noise be 
small perturbations. 

Before proceeding with the formal derivation, let us first point out some subtleties related to the stochastic case 
as opposed to the deterministic one. In the latter, it is customary to define a curvilinear coordinate system (s,r) in 
which r = stands for the curve representing the front position, which can be associated for instance to a level curve 
of the appropriate field. The scheme assumes that the front thickness is small compared to the radius of curvature, 
and that the relaxation of the internal degrees of freedom of the front is fast compared to the time scale of the 
long wavelength deformations of the front. When noise is present in the field equation, the appropriate curvilinear 
coordinate system cannot be associated to level curves, which are very rough at length scales smaller than the front 
thickness. On the other hand, at larger scales, which are the ones we are interested in, a coarse grained description 
makes perfect sense. This is actually implicit in the very idea of the stochastic eikonal equation. One can think 
of different schemes to explicitly define such coarse-grained description, all of them equivalent. However, since the 
rest of the derivation cannot be carried out explicitly in full rigor in any of those, and since the result is expected 
to be independent of the details of the definition of the coarse-graning, we will proceed more or less heuristically. 
In essence this is a reformulationn of the approach introduced for the derivation of the diffusive wandering of fronts 
in one dimension discussed in Ref. E^ j. There, the basic idea was that only the low-frequency components of the 
noise are responsible for the front wandering so the high-frequency components can be implicitly integrated out. The 
effect of the high-frequency components of the noise is thus to renormalize the mean front profile. As a consequence, 
they renormalize the front velocity, and also the diffusion coefficient. More precisely, this means in our case that the 
high-frequency renormalization is carried out by the Stratonovich term in the function h. Once this term is explictly 
extracted, the remainder is an Ito multiplicative noise. Then the high-frequency components of this term are irrelevant 
and can be averaged out, while the low- frequency components will be responsible for the roughnening of the front at 
scales larger than the coarse-graining length. In the case of Id fronts, the resulting diffusion coefficient for the front 
wandering has been rederived rigorously in Ref. [^5| . Unfortunately, the identification of the collective variable in Id 
cannot be immediately generalized to higher dimensions, so we must still rely on a less precise formulation and check 
the consistency with rigorous results and numerical simulations a posteriori. 

After the above considerations, we can make the following theoretical construction. We assume we have solved 
the field equations with noise without any approximation. We now coarse grain the fields with some local average 
procedure, both in time and space, and use the coarse-grained fields to define a curvilinear coordinate system based, 
for instance, in terms of level curves at any time. At short length and time scales this coordinate system is smooth 
and, in principle, we could write the full field equation (still with the bare fields) in these curvilinear coordinates. In 
the absence of noise, an expansion in the front thickness would unambiguosly yield the terms which are dominant 
in the range where the eikonal equation is devised for, namely for sufficiently long length scales (small curvatures) 
and long time scales. Terms such as second derivatives on s and the time derivative would automatically drop out 
of the description. In the presence of noise stochastic case, this is not so automatic unless the fields themselves are 
coarse-grained so that they are also sufficiently smooth in space and time. We assume that, for the coarse-grained 
fields, the order of the different terms in the front thickness will be the same than for the deterministic case, when such 
expansion makes sense (excluding pulled fronts, for instance). We claim that such assumption is the one implicit in 
the very idea of the existence of a stochastic kinematic formulation of the front dynamics. Then, for the coarse-grained 
fields, Eq. (p3) is expected to reduce, in analogy to the deterministic case, to 
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+ £>K(s,t)^ +h(cf>) + v n (s,t)^ + e^ 2 n(cj);r,s;t) = 0, (9) 

where n is the local curvature and v n the normal velocity of the front. This normal velocity provides the evolution 
of the curvilinear coordinates in which Eq. (j5|) takes the form of Eq. (^|), and is the fundamental quantity we 
are interested in. The noise term in Eq. (^|) must also be considered as coarse-grained, with the high-frequency 
high-wavenumber components integrated out. 

At this point of the derivation it is useful to consider the Id problem which corresponds to the zeroth order of the 
eikonal description. This is defined by neglecting the curvature and the fluctuating term in Eq. (^): 

° = ^+^A)^+h(0 o ). (10) 

This is the eigenvalue problem that gives the renormalized velocity v in the Id problem as obtained in Refs. [ ]l2| , |l3| |. 
Note that this equation does contain noise effects through the high-frequency renormalization provided by the 
Stratonovich term. In fact, the effective velocity v(e\) resulting from Eq. (|l^) has contributions from all orders 
in the dressed noise intensity e\. An explicit first order in E\ approximation will be given in the Appendix. In the 
problem defined by Eq. (|^) curvature and fluctuations will be taken as small perturbations. Hence we assume for the 
field <p and the velocity v n expansions of the form 

4>(r,s,t) = 4> (r)+6cjy(r,s,t), (11) 



v n (s,t) = v(e\) + {3(s\)k(s, i) + Sv(s,t), 



(12) 



where <fia(r) and v(e\) are the solution of the Id problem of Eq. jltj|). The term /3«;(s,f) is a curvature correction 
and Sv(s,t) describes fluctuations. Linearizing in both perturbations S(p(r,s) then verifies 



where 



= f 8<f> + (/3k + 6v)^p- + Dn^p- + e 1/2 ft(0 o ; r, s; t), 
Or Or 



V = D d^ + ^ ) dr + d4> 



0=00 

Taking the derivative of Eq. (|l^) with respect to r, it is a simple matter to prove that 

d(p 



u = 



Or 



(13) 



(14) 



(15) 



is the right eigenvector of T with zero eigenvalue. Due to the non-hermiticity of T, finding an analytic expression for 
the left eigenvector, u°, is not trivial. Notice that because of the vectorial character of the field, the simple expression 
obtained in Ref. |13| for scalar fields is in general not applicable. Nevertheless, the corresponding eigenvector can 
always be obtained at least numerically. 

Now, taking Eq. ( p"3] ) and performing the scalar product with u°, we obtain 



k(u°, Duo) + 0k(u°, u ) + e 1/2 (u°, ft(0 o ; r, s; *)) + 6v(s, t)(u°, u ) = 
where the scalar product is defined by 

(f,g)=£ J drf t {r) 9l {r) 

i 

Owing to the independence of the two first-order perturbations (curvature and fluctuations), we get 

(u°,Cu ) 



P(e\) = - 



(u°,u ) 



(16) 



(17) 



(18) 



and 
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MM) = - £ 1/2 ^|^. (19) 

The stochastic process ( |l9| ) is not white since the high-frequency components of f2 have been integrated out by 
the coarse-graining procedure. However, we can now restore them harmlessly by replacing H with the original 
multiplicative white noise process. By doing so, we are modifying the part which is not intended to be accounted for 
by the very eikonal description, while the treatment is simpler. Analogously, once the dependence on the cut-off A 
has been explicitely worked out, we can take the process (|l^) as delta-correlated in space. We will explicitly check the 
limit of validity of the eikonal description here proposed when scales comparable to the front thickness are reached, 
in the sections below. 

The resulting stochastic eikonal equation with the explicit dependence on the original noise parameters then takes 
the form 

v n (s, t) = v(e x ) + P(s x )k(s, t) + D)' 2 {e, e x )C(s, t) (20) 

where v(e\) is defined by Eq.(|l0|), /3(e\) is given by Eq.(|J), and the noise £(s, t) is a zero mean Gaussian white 
process with 

(as,t)as',t'))=26(s-s')6(t-t% (21) 
which follows from the statistical properties of f2 with 

Df(e iex )=e (u O iUo)2 ■ (22) 

Note that the dependence of Df on the dressed noise intensity s\ comes from the dependence on the same quantity 
of (po (and hence of uo and u°) as solution of the renormalized problem of Eq. (fio|). 

The above equations constitute the main result of the first part of the paper. Although the derivation is not 
rigorous because the coarse-graining could not be carried out explicitly, the result is appealing from the theoretical 
point of view in that it separates the problem into an effective deterministic one, where the original field equations are 
modified by additional deterministic terms which depend on noise parameters, plus an additive noise which would be 
the necessary one to describe the wandering of the problem in one dimension. For the renormalized one-dimensional 
deterministic problem, therefore, the two present perturbations, curvature and noise, decouple from each other. 

An important point to emphasize here is the separate dependence of the result on two noise parameters, namely s 
and £\. While the renormalized velocity and the coefficient of the curvature depend solely on ex, the effective noise 
intensity D f depends separately on both. This clearly illustrates how the ultraviolet cutoff is an additional parameter 
of the problem when the noise is multiplicative, in correspondence to the fact that the continuum limit is not well 
defined for noise delta-correlated in space. It is also important to remark that our derivation procedure is expected 
to be valid for small noise intensity e, but contains all orders in ex- This was already the case in the one-dimensional 
derivation of Ref. fll3fl , where the small noise approximation was phrased in terms of a separation of time scales. 
The connection between that scale separation, the coarse-graining procedure and the small noise expansion has been 
recently clarified in Ref. |Q , where a rigorous derivation of the result of Ref. p3] has been presented for the case of 
a single-component field in one dimension, in terms of suitable projection techniques. Unfortunately that rigorous 
derivation is based on the identification of a specific collective coordinate which has no simple extension to higher 
dimensions. Nevertheless, the fact that this approximate procedure has proven correct in Id gives further support 
to our main result above, which is not claimed to be rigorously proven. In the following sections we will check 
this prediction against analytical results and numerical simulations of the full reaction-diffusion equations in explicit 
examples. We will see that, the dependence on the cutoff A is essential not only for a quantitative description of the 
problem, but is crucial to predict nontrivial phenomena such the transition to pulled fronts, in which the whole eikonal 
description fails. This failure of the present description is signaled by the vanishing of the effective noise intensity Df 
when that point is reached. In fact, Df is linear in e to lowest order, but has a complicated dependence on ex- As we 
will explicitly see, the partial resummation of orders in ex captures important physical features of the problem. For 
instance, it allows the non-monotonic dependence and eventual vanishing of the front diffusion coefficient Df. We 
expect the pushed-pulled transition to occur exactly at this point. Another qualitative change captured by the above 
resummation is the destruction of the front itself, associated to the fact that the front thickness may diverge in some 
circumstances. This phenomenon is the signaled by a dicergence of Df at some finite value of ex- 

Although the predictions above are expected to be correct for ex ~ 1 as long as e <C 1, in practice this may be 
limited by the fact that v(ex) and the Goldstone modes are not in general analytically known. In such case one can 
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rely on numerical resolutions of the eigenvalue problem posed by Eq. (|o|), or alternatively one can find v(e\) as a 
systematic expansion in powers of E\ as described in the Appendix. On the other hand, it is worth remarking that (3 
can be a nontrivial function of e\ only for multicomponent fields, that is for pulses. For one-component fields (fronts) 
the coefficient f3 is not renormalized by noise. 

We remark that this derivation is meaningful only when pulses are involved or, in the case of fronts, when the 
relevant dynamical regime is the pushed one. In the last Section of this paper, we shall give more details about 
the main differences between pushed and pulled fronts, and about their different response to noise. For the time 
being, we point out that the lack of time-scale separation between the relaxation of the zero mode and the other 
eigenmodes of the spectral operator for pulled fronts prevents in general from constructing a local equation for the 
interface motion ]29| ] . Examples of reaction-diffusion systems which do not admit local kinematic descriptions are the 
phase-field formulations of solidification |3(]j4(| or viscous fingering p2l. In those cases a nonlocal interface equation 
does exist so an extension of our derivation is in principle feasible. This may be particularly interesting in cases such 
as in Refs. |40|]32[ | were the relevant fluctuations may be external. 



III. KINETIC ROUGHENING AND CONNECTION TO EQUILIBRIUM FLUCTUATIONS 



In the previous Section, we have derived a stochastic sharp interface approximation for a generic RD system with 
multiplicative noise. A pictorial description of a noisy front is given in Fig. u\ 




FIG. 1. An example of a noisy front with a level curve which defines the precise location of the front. 

We see in this figure that noise in the RD system induces fluctuations in the front shape, thus generating roughening 
of the sharp interface that should emerge at the eikonal level of description. The identification of universality classes 
of kinetic roughening will come naturally at this level of description. In particular we will establish the connection 
with the universality classes defined by the Kardar-Parisi-Zhang (KPZ) Equation |Q and by the Edwards- Wilkinson 
(EW) Equation @. 

The stochastic eikonal equation ( p(j| ) is written in instrinsic, rotation invariant form. For the purposes of scaling 
theory it is convenient to write it in cartesian coordinates. The front location is then given as y = h(x, t). Retaining 
only the relevant nonlinear terms in the Renormalization Group (RG) sense, we then recover the KPZ equation 

Bh d 2 h A fdhV , . 

with 

(H(x, t)v{x',t')) = 2D KPZ 5(x - x')S{t - t'). (24) 
The KPZ parameters turn out to be related to the eikonal ones as 

D KPZ = D f , v = p t X = v. (25) 
In the special case of v = 0, the EW equation is obtained, 
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dh d 2 h 

with 



St -^ + K*,*), (26) 



{li(x, t)fi(x', t')) = 2D EW S{x - x')S{t - t'), (27) 

and again 

D EW = D f , v = p. (28) 

These equations are well known to be the paradigm for many different growth processes (]33| , j34| . Even if the microscopic 
dynamics of the system under study may correspond to different equations of motion for the respective interfaces or 
surfaces, nevertheless the KPZ and EW Equations do capture the universal features of the system, namely the scaling 
properties of fluctuations. 

Usually, such effective equations cannot be derived from the original microscopic description of the particular 
systems and are introduced on a phenomenological basis, relying on the claim of universality within a RG framework. 
Nonuniversal quantities such as prefactors of scaling functions, affected for instance by the noise intensity in the 
interface equation cannot be derived. In our case, we are able to compute the noise intensity in the eikonal equation 
so we can also predict the nonuniversal prefactors if the noise is known in the reaction-diffusion level of description. 
For instance, not only the scaling of the interface roughness with system size can be predicted, but also the actual 
values of average interface roughness in terms of the original microscopic parameters of the RD model are worked out. 

As a test of our derivation we will now check consistency with equilibrium fluctuation theory. The connection be- 
tween bulk thermal fluctuations and fluctuations of the interface between thermodynamical phases can be established 
rigorously in the case of equilibrium fluctuations. This is possible because a free energy functional does exist and the 
sharp-interface limit can be performed at the level of the free energy itself. Then, the fluctuations can be obtained 
independently from the free energy at each level of description (either bulk or interface fluctuations), consistently with 
the fluctuation-dissipation theorem. The important difference here is that no dynamical equation must be invoked 
but only equilibrium properties. On the contrary, in the more general case where there is no fluctuation-dissipation 
theorem and not even a free energy functional, we must rely on dynamical equations. In the case of equilibrium 
fluctuations, however, we must reproduce the known correct result. As we will see below, this case falls in the EW 
universality class. 

The general solution of the EW Equation, Eq. (p6|) , is known H]. Consider the interface h(x,t) and its discrete 
Fourier Transform h q (t), defined through 

h(x, t) = h(t) exp(igx). (29) 
i 

It is possible to show J34| that the long time limit of the spectrum S(q,t) = (\h q (t)\ 2 ) is given by the expression 

lim S(q,t) = j ^f. (30) 



Our strategy now will be to calculate explicitly the spectrum (|3fj) in terms of the coefficients predicted for our eikonal 
equation, and then show that the resulting expression coincides with the independent result that can be obtained 
from equilibrium fluctuation theory. Hence we insert an additive noise in the original RD system, that is in setting 
.9 = 1: 

^=OT 2 ^ + ^) + ^(x,f), (31) 

with 

(r } (x,t) V (x',t'))=26{x-x , )6(t-t'). (32) 

We consider a F((j)) with a symmetric double- well form, i.e. the deterministic part of Eq. (jjH]) has a Id solution with 
zero velocity cj> = <po(x). This is what corresponds to the usual time-dependent Gizburg-Landau Langevin equation 
for a nonconserved order parameter (Model A in the Hohenberg-Halperin classification pq|), where noise intensity 
must be identified as e — ksT. Since g = 1, our expression for the noise intensity at the sharp interface level takes 
the simpler form 
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Df = j, (33) 

, ( d<pa 
ax 



dx 

which, according to (|30|) and performing the identifications of Eq. fl28| ) with (3 = D, produces 

Hm (\h q (t)\ 2 ) = ' - n2 ^- (34) 

LDj^dx^t) q 

On the other hand, we can take the sharp-interface limit on the free Ginzburg-Landau free energy. The calculation is 
standard (see for instance Ref. |37f]) and yields the interface free energy 




Fi=a dx\ 1+ — , (35) 



where the parameter a is identified as the intcrfacial tension and can be evaluated from the bulk free energy of the 
system |$7| as 

(36) 

where 0o is the corresponding kink solution. For soft (long wavelength) deformations of the interface, the excess free 
energy reads 

2 



Oh 



AF ^2j dx {8x)> < 37 > 

and the corresponding stationary spectrum of fluctuations, consistent with the fluctuation-dissipation theorem takes 
the form 37 

(Kl 2 ) = (38) 

Using Eq.(|36|), this expression yields the same result of Eq.([34j). This proves that the front roughening obtained 
from our derivation is exact in the case of equilibrium fluctuations, and by extension in the additive noise case. We 
have shown this for the case of a non-conserved order parameter. In the conserved case (Model B of Ref. pf|), the 
projection to a sharp interface description yields a non-local equation and therefore lies outside the validity of our 
theory. Similarly, equilibrium fluctuations have also been studied in both sharp-interface [B8| and phase-field j39 



formulation in the context of solidification, which also yields non-local interface dynamics. The universality classes in 
those cases are not well established. 



IV. APPLICATION TO A PROTOTYPE MODEL OF FRONT PROPAGATION 

To illustrate our general theory we consider here as an example the propagation of a scalar front 

^ = DV 2 + F(0,a)+£ 1/2 <?(0>7, (39) 
with the noise correlator defined as in Eq. (^|). We specify our prototype model through the following definitions: 

F(0)=0(l-0)(0 + a) (40) 
g{4>) =0(1-0), (41) 

and we will consider a front of the 0=1 state invading the = one. The constant a is a control parameter. As it 
is well known for fronts without fluctuations, the deterministic force given by ( f40| ) leads to different modes of front 
propagation depending on the value of a 0. The front velocity depends also on a. The choice of the coupling function 
( pl| ) for the multiplicative noise term of Eq. (^) is the simplest one that preserves the stationary states = and 
0=1. The fact that the noise term vanishes in the two asymptotic states prevents nucleation phenomena in the 
invaded state. The form of g(4>) is also such that the multiplicative noise term arises naturally as the fluctuation of 
the control parameter a. Moreover, for the prototype model here proposed, the corresponding function /i(0) defined 
in Eq.(^) which appears in the renormalized equation takes the same functional form as /(0), only with renormalized 
coefficients. This allows a simpler analytical treatment and intepretation of the results. 
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A. The Id case revisited 



This prototype model has already been analyzed for Id in Refs. [|r2j,|l3j], where in the regime — e\ < a < 1/2 — e\ it 
is proven to result in a renormalized average velocity 



2a + 1 



and in the diffusion coefficient 



D * 



n^e 2 ^ (#oM)V(0o) 



(42) 



(43) 



Here </>o is the solution of the zeroth order equation: 



__ + ,_ + ^ o) = . 



(44) 



Notice the dependence of both these quantities on A through the s\ parameter. Specifically about the diffusion 
coefficient, the first order dependence in e is apparent, while the functions present in the integral and defined through 
(fl4|) contain all orders in e\. This is in contrast to the renormalized velocity, which depends solely on ex. This means 
that, both s and A can be determined independently from separate measurements of both the ballistic and the diffusive 
components of the front propagation. This is quite remarkable since it provides indirect means of measurement of the 
(microscopic) noise, which may not be directly accessible in many cases. 

As an illustrative and well-controlled example, we have explicitly tested the prediction of the dependence on A of 
both velocity and diffusion coefficient with direct numerical simulation of the RD equation. The first set of results 
is shown in Fig. |^. Here data from a simulation of the RD equation with the reaction term ( pp[ ) are reported and 
compared with the theoretical prediction (^2|) for two different values of A = 1 and A = 5, in some dimensionless units 
of the simulation. 




1 = 5 



0,5- 



°0 ' 50 ' 100 150 200 

t 

FIG. 2. Change with A of the average renormalized velocity. Values of the parameters are: a = 0.1, e = 0.1, A=l,5. The 
theoretical values of the corresponding average velocities as from Eq. (^), 0.948 and 0.866, are also plotted. 

Notice that the result is indeed sensitive to the microscopic cutoff A of the noise, even though this length is 
significantly smaller than the front thickness, of the order of 25, in the same units. 

Regarding the diffusion coefficient (^), we measured the mean square displacement of the front position. If we 
define the front position as 

/•OC 

z(t) = / dx<f>{x,t), (45) 

J xo 

then the diffusion coefficient (^) is related to the mean square displacement 
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V(* 2 ) w 



(46) 



as 



A 2 {t) ~ 2D ft 



(47) 



In Fig. 



the quantity A(i) is plotted 

10 




i 



FIG. 3. Change with A of the diffusion coefficient for the front wandering in Id. Parameters have the same values as in Fig. 



while the values of A run from 1 to 15. The diffusive 



The values of the parameters are the same as in Fig. 
behavior and the dependence on the value of A is manifest. 

In this particular case of Id, a more systematic derivation of the diffusion coefficient has been reported in Ref. [fL5| . 
By proper identification of the natural collective variable which describes the front wandering as strictly (not just 
asymptotically) diffusive, it has been shown that, in fact, the result first found in Ref. |l3| is rigorous to first order in 
e. Together with the case of equilibrium fluctuations, this is the second rigorous test of our general theory. 



B. The 2d case 

Let us now consider the propagation of a front in the prototype model of Eqs. (^, |{], |l]) in 2d. We have already 
shown that in two dimensions the eikonal equation reduces to either the EW equation or the KPZ equation, depending 
on whether the velocity of a planar front is zero or non-zero respectively. In order to make a numerical check of our 
theory we consider then the simplest case of zero velocity, with expected EW scaling. This case corresponds to the 
choice a = —1/2 in Eq. (f40|). As seen explicitly in Eq.(|4^) it turns out that the renormalized velocity is also zero, 
since the noise does not break the symmetry associated to the double-well form of the deterministic potential. It 
is important to remark that, unlike the Ginzburg-Landau model discussed in Section III for the equilibrium case, 
the noise is now multiplicative, and no fluctuation-dissipation theorem can be invoked. Hence the first-principles 
derivation of the fluctuation spectrum in the sharp-interface description is no longer available. 

We thus rely on the dynamical equation 

^ = DV 2 4> + 0(1 - m - i) + eVVtl - 0Mx, t), (48) 

where the noise r\ is defined through the usual correlator of Eq. (0) . 

Now, our basic goal is to connect this level of description with the eikonal level, determining thereby the noise 
corrections to this equation. Therefore we assume a stochastic eikonal equation of the form 

/ Vn ( S ) = -p K+ D i / 2 a Sl t) ( 

\({(s,t){(s>,t>))=2S(s-s>)5(t-t') 1 ' 

where the coefficients (3 and Df are given by 
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Notice that we have directly taken into account that the bare as well as the noise-renormalized velocities are both 
zero (see for example Also notice that the renormalization of the curvature term is absent for both e and e\ 

expansions due to the fact that we are now dealing with a front (N = 1). 

Now, in order to calculate Df, we need to specify the solution of the Id model. This is known and for a = —1/2 is 
given by 

</> = ^(1 -taahkx), k = (51) 



and 



u vo ,,„>. 

u=u = —. (52) 



The integral in (pOj) can then be computed exactly and gives: 




-35yi-2e/A2- (53) 

The above result clearly illustrates the different treatment of the parameters e and e/X 2 . The result is first order in 
e and contains all orders in e/X 2 . It is interesting to remark that the partial resummation of all orders in e captures 
important nonpertubative phenomena. For instance, the divergence of Df at e/X 2 = 1/2, reflects the fact that at this 
point the front itself is destroyed, or, equivalently, the front thickness becomes infinite. This is equivalent to reaching 
a critical point, except that this is not the equilibrium one because the noise is multiplicative. Remarkably, our result 
for additive noise case does not capture that feature because it is only first order in the noise strength. In the case 
with an asymmetric double-well potential, which has a finite front velocity, it was explicitly checked numerically in 
Ref. Jl3| that the diffusion coefficient of the front has a non-monotonic dependence with e/X 2 . Most importantly it 
vanishes at a finite value of e/X 2 which corresponds exactly to the point were the front reaches the pushed-pulled 
transition. Again we see that the resummation of orders e/X 2 captures important physical information (see discussion 
in Section V). 

We now come back to the numerical test of the EW scaling of our particular symmetric case, with the identification 
(p8|). Accordingly, we can rewrite the complete power spectrum as it is known theoretically [pi| in terms of the 
coefficients that we just calculated, and compare it with data from a direct simulation of the field model Eq. (fl8|) . 

From |34j the power spectrum S reads, 

S(q, *) = (1 - exp(-2*Vt)). (54) 

L vq z 



In terms of the parameters in the original field equation this yields 



1 D 



S(q,t) = j-pL(l-exp(-2Dq 2 t)), (55) 

where Df is given by Eq. (|53|) . Accordingly, the fluctuations of the front position in the RD model Eq.(|39|) for length 
scales larger than the front thickness itself must obbey the spectrum defined by Eq.(|55|). It is worth remarking that 
the prediction is not only for the universal features, namely the shape of the scaling function and the exponents, but 
for the actual absolute values of the spectrum. We have performed simulations of the field RD model with a correlator 
of the form 

fo(x, t)r,(x',t')) = ^(1 - |x- x'|/A)(l - \y - jfl/A)0(l - |x - x'|/A)6(l - \y - y'\/X)6(t - t') (56) 

This corresponds to assuming that at any time, the noise takes the same value in a square of side A, uncorrelated 
with the neighboring squares. This is done for simplicity but no significant dependence is expected on the details of 
the shape of the spatial correlation provided A is kept smaller than front thickness. 

We have studied the fluctuations of the internal level curve of the front <f> — 1/2. In Fig. |] we show the scaling 
region, with the correct slope and location of the curve. More remarkably, Fig. [s] shows the measured spectra for the 
simulation of the field equations compared to the prediction given by Eq. (B5J) . 
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FIG. 4. Numerical data from 2d simulation of the RD field equations for the protorype model with a — —1/2 and analytical 
prediction for the power spectrum (^34|) in the scaling region. The time is here 1000, and the parameters of the simulation are 
L = 100, e = 5, A = 5. 




FIG. 5. Numerical data from 2d simulation of the RD field equations for the protorype model with a = —1/2 and analytical 
prediction for the power spectrum (^). The three sets of data points refer to times 50, 250 and 500. The parameters of the 
simulation are L — 500, e = 1, A = 2. The value of q c corresponds to a wavelength of the order of the front thickness. The 
analytical prediction is only intended for q < q c . 



It should be stressed that in this comparison there is no free parameter. It is also interesting to observe the 
deviations from the prediction at length scales smaller than the front thickness. In this high-q region the data also 
collapse but not to EW scaling. An estimate of the exponent a in this region is around 3/2. 

For values of a such that the front has a finite velocity, one would expect that the scaling would be given by that 
of the KPZ universality class. Although the scaling function in that case is not exactly known, the prediction of our 
eikonal equation is expected to fit the data for the corresponding RD model also without free parameters. We have 
not checked this case because it is obviously more involved and less conclusive because of the practical difficulties to 
reach the scaling regime already at the eikonal level of description |4(J| . 

From a practical point of view it is to be remarked that the noise intensity cannot be increased arbitrarily in a 
simulation without destroying the front itself. This can be easily seen in a numerical simulation. Although the noise 
vanishes asymptotically in the stationary states <p = and cj> = 1, if noise is sufficiently strong it may be capable to 
nucleate the other state in the region not too far from the front. We have found that this effect is more pronounced 
in the region behind the front <p = 1. For a given A there will typically be a maximum value of e after which the front 
is essentially destroyed. If we increase A the effect is milder. In order to make the front roughening appreciable in 
not too large system sizes, it is thus convenient to have a moderately large A, which in turn will allow larger values of 
e. Typical values that we have considered are in the range of A = 2, 5 n units for which the front thickness is of the 
order of 25. 
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V. THE PUSHED - PULLED TRANSITION 



As we have being mentioning so far, two classes of fronts must be distinguished from a dynamical point of view, 
the so-called 'pushed' and 'pulled' fronts f|,^j| . The simplest is pushed case, in which the front propagation depends 
on the full non-linear structure of the equation of motion and the front is said to be 'pushed' by its internal part. 
This is usually the case when the invaded state is locally stable. On the contrary, if the invaded state is unstable, it 
can happen that the relevant dynamics takes place in the semi-infinite leading edge region ahead of the front itself. 
Then the propagation of the front is governed by the growth and spreading of linear perturbations in that region 
which 'pulls' the front. In this case the linearization about the unstable state accounts for its dynamical behavior [Q, 
but there is degeneracy of propagating velocities @]. In the present context, the most important distinction between 
the two situations is that, while for pushed fronts, the relaxation of bulk modes is exponential, for pulled fronts it 
is algebraic, as a result from the fact that the linearized operator describing perturbations around the stationary 
propagating mode is gapless. This means that for pulled fronts, there is no natural time-scale separation which allows 
for the decoupling of the interface modes from the bulk ones. Our derivation, and the whole idea of a kinematic 
moving-boundary approximation, is not valid for pulled fronts |29| ]. Nevertheless, we will show that our theory does 
predict correctly its failure at the pushed-pulled transition. 

The intrinsic differences between pushed and pulled fronts also show dramatically in the statistics of fluctuations. 
In Id, for instance the front wandering turns from difussive (pushed fronts) to subdiffusive (pulled fronts) |l3|]l5| . In 
turn, pulled fronts in 2d with multiplicative noise have been found to belong to a different universality class than the 
ordinary KPZ |]l6| , |l7| , as opposed to the KPZ scaling of pushed fronts. 

In Ref. Jl^] it was already observed that the diffusion coefficient of fronts in Id was a nonmonotonic function of e/A 
which vanished at a finite value of that parameter. After that point, subdiffusive behavior was found. The picture was 
confirmed and completed in Ref. ||l5f . Here we want to stress that the point where that transition occurs corresponds 
exactly to the pushed-pulled transition. In fact, the renormalized equation defined by Eq.(^) for the protoype model 
takes the same form as the original one, with renormalized coefficients. The deterministic equation has the transition 
to the pulled regime at a = 1/2. The same transition in the renormalized equation does occur at a = 1/2 where 
a = a + s\. It is thus clear that, for parameter values of the deterministic equation in the pushed regime, increasing 
noise intensity e or decreasing A will imply a transition to the pulled regime. While the front velocity will not 
experience any dramatic effect, the fluctuations (wandering in Id or roughening in 2d) will be dramatically affected, 
since the universality class will change. In Fig. [| we show the change from diffusive to subdiffusive behavior induced 
solely by a change in the effective noise instensity e\. 



— £ = 0.5 

£ = 0.1 




\o 100 

t 

FIG. 6. The pushed-pulled transition in Id. The system size is L = 2000 and averages have been carried out on 3600 
realizations of noise in the pulled case and about 1000 in the pushed one. 

Once more it is remarkable the dynamical importance of the noise correlation length A in the long-wavelength 
behavior of the front. The same effect should be expected in 2d, namely, the scaling could be KPZ or non-KPZ 
depending solely on noise parameters. To our knowledge, this is the first time that such a dramatic effect of noise is 
reported. Remarkably enough, while our eikonal description is not able to describe the second regime, it does predict 
the transition at the right values of parameters. 
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VI. DISCUSSION AND CONCLUSIONS 



The formulation in terms of kinematic cikonal-like equations provides a useful framework for studying front or pulse 
propagation when one is interested in long spatial and temporal scales. This kind of equations have advantages both for 
numerical simulations and for theoretical analysis, and has fruitfully been used for instance in the iphcnomcnological 
modeling of excitable wave propagation in disordered and noisy media, by the ad Ziocprocedure of adding fluctuations 
to a generic eikonal equation. In this paper we have derived stochastic eikonal equations from the more microscopic 
RD field equations with noise, which can be multiplicative. 

The derivation presented here relies on the hypothesis of separation of scales between the front dynamics at large 
scales and the internal degrees of freedom of the front. That means that it is not valid for pulled fronts, which are 
indeed known not to have a local, eikonal-like description even in the absence of noise. However, even for the pushed 
case, the usual projection techniques for derivation of sharp-interface equations cannot be simply extended to the 
stochastic case due to the fact that the noise contains the full range of length and time scales. Hence our derviation is 
formulated within a coarse-graining scheme, which we claim makes sense for the kind of problems which admit a local, 
cikonal-like equation. Derivation of sharp-interface approximations with noise have been possible so far for thermal 
noise only, in systems with local equilibrium. For multiplicative, generically external noise, however, the absence of 
a free energy and a fluctuation-dissipation theorem requires an alternative scheme based on dynamical equations. 
We have explicitly checked that our scheme is exact for the cases of equilibrium fluctuations and also reproduces a 
rigorous result for multiplicative noise in d. For the general case, however, we rely on numerical tests to fully justify 
its validity. In any case, the excellent agreement with numerical and analytical tests clearly suggests that a more 
rigorous derivation should be possible. An extension of our procedure is also conceivable in RD systems for which a 
nonlocal sharp-interface description exists for instance in solidification or viscous fingering. 

One of the main points of this paper has been to clarify the role of the spatial cutoff of noise correlations A. On the 
basis of recent rigorous mathematical findings on the spatio-temporal white-noise limit, we have argued and shown 
in explicit examples, that the common splitting of the Stratonovich white noise in a term which is singular as A — > 
plus an ltd noise, does capture the correct dependence on A of macroscopic quantities (at scales much larger than A). 
The remaining spatially-correlated ltd noise will only carry a weak dependence on A, noticeable at length scales of 
the order of A, and can be treated perturbatively. The A — > limit can thus be taken safely for the remaining Ito 
noise, once the singular part has been extracted. Within this scheme, macroscopic quantities such as front velocity 
or front roughening properties, depend separately of the noise strength e and on e/X d . Incidentally, this implies the 
possibility of measuring the (microscopic) noise parameters e and A from the macroscopic dynamics of the fronts. 
Most remarkably, a partial resummation scheme can be naturally defined for the parameter e/X d which captures 
some important, nonperturbative physical phenomena which would be missed otherwise, such as the transition from 
pushed to pulles regimes, or the destruction of the front associated to a divergence of the front thickness. These 
phenomena, at the edge of validity of the sharp-interface approximation, are detected respectively by the vanishing 
and the divergence of the effective noise-strength on the eikonal equation. 

Eikonal stochastic equations like the ones derived here are directly related to the EW and KPZ equations of kinetic 
roughening. By using known results for these equations and the results presented here one can predict roughening 
properties of noisy pulses or fronts appearing in field equations. We have used such correspondence to check the 
predictions in a case of zero velocity and additive noise, for which the general theory of equilibrium fluctuations 
can be directly applied to the RD equation. Analytical results obtained for the corresponding (Edwards- Wilkinson) 
eikonal equation have been identical (including prefactors) to that independent calculation for the RD system. 

We have also applied our results to a prototype model with multiplicative noise constructed to represent a variety of 
different front propagation regimes by changing a single control parameter. While this model has already been used to 
study effects of noise on Id fronts, we have addressed here a more specific 2d effect such as the case of front roughening. 
In a zero-velocity case, simulations of the reaction-diffusion equation have presented a very good agreement with the 
predictions of roughening for the corresponding EW equation with no adjustable parameter. The results show that 
the dependence on A is quantitatively important even if A is significantly smaller than the front thickness, which 
may seem counterintuitive. Although we have not checked the analogous results for a nonzero velocity, which would 
correspond to the KPZ equation, we expect our results to be valid also in this case. 

Most interestingly we have explicitly checked the prediction of qualitative changes as e/X d is varied, such as the 
transition from pushed to pulled propagation regimes. We have directly observed the change in the wandering exponent 
in Id as A is decreased. The immediate extension of this result implies that a transition from KPZ to non-KPZ scaling 
is to be expected in higher dimensions. We thus conclude that the dependence on the spatial cutoff of the noise may 
have dramatic effects, not only on the nonuniversal quantities but also on the universal ones. 
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In summary, we have derived stochastic eikonal equations from stochastic RD equations completely specifying the 
parameters of the noise. Although this is not a systematic derivation, we have presented a wealth of independent 
evidence on the validity of our results both from analytical calculations and from numerical simulations of different 
systems and looking at distinct noise effects, both qualitative and quantitative. These results should be of interest 
for theoretical and practical purposes both in the study of kinetic roughening and in the context of propagation of 
chemical waves. 
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APPENDIX: SYSTEMATIC EXPANSION 

In some circumstances, the eigenvalue problem defined by Eq. (|n]) may not be solvable analytically, due to the 
functional form of h(<po), which has been renormalized by noise, even though the 'bare' problem may be solvable. In 
such cases it may be useful to solve for the case e = (no noise) and obtain corrections in a systematic expansion in 
e. This is the aim of this appendix. 

Let us consider then the zero noise eigenvalue problem defined by the (deterministic) Id equation 

~d 2 (j>T> d<f)T> 

= D + V Q^ hf(0x>) (57) 

or^ or 



We now expand the fields of the noisy 2d problem of Eq. (J9j) as perturbations in e and k of the solution of Eq. (J57 

M) = <Mr) + <ty>(r, M) (58) 



v n (s, t) — vq + as + f3n(s, t) + Sv(s, t) 



(59) 



We get at the linear order 



= f 8<j> + eG{<hj) + (as + f3k + 5v)^- + t)k^- + e 1/2 tt(<^; r, s; t) 

or or 



where 



. ~d 2 d df 
or z or ocp 



<t>=4>- 



(60) 



(61) 



x> 



Now the right eigenvector of f reads 



d(f>T> 
U (0) = ~gy~ 



(62) 



which is independent of both e and k. Taking (|o|) and performing the scalar product with the left eigenvector u^°\ 
we get 



fc(u<°\ .Du ( o)) + (as + /3fc)(u(°\ u (0) ) + £ ( U ( 0) , G(«fr>)) 

V°\ n(cf>v; r, s; t)) + 5v(s, t)(u (0 \u {0) ) = 



(63) 



As curvature and fluctuations are linear perturbations we get 
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(u(°),G(<M) (n (0) ,Ou ( o)) 

(u(°),u (0 )) ' ; (uW,u (0) ) 

and 

,(o) 



c , ,x 1/2 (uW,n(^; M ;t)) 

dw(s,i) = -e 7 7 fTvi \ • 

(u (0) ,U(o)) 



The stochastic eikonal equation is now 

v n (s, t) — v + ae + /3k(s, t) + D 1 f /2 (e)C(s 1 1) 
where the noise ((s, t) is a zero mean Gaussian white process with 

(C(s,t)as',t , ))=25(s-s')6(t-t% 

and 



, , /drX)ij«i O) «j(o)fli(0i>)fc-(^ 
Dfi£) = £ (u(o),u (0))2 
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